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I .  SUMMARY 

This  report  presents  a  study  o£  the  nature  of  the  unsteady 
forces  in  the  field  of  a  propeller  rotating  In  the  vicinity  of  an 
appendage.  The  propeller  Is  assumed  to  be  one  of  a  high  aspect  ratio 
while  the  appendage  although  of  finite  width  Is  assumed  to  be  Infinitely 
long  and  thus  the  problem  Is  reduced  to  a  study  of  the  two-dimensional 
flow  field  around  two  flat  plates.  An  essential  feature  of  the  analysis 
Is  that  the  mutual  Interference  effects  of  propeller  blade  and  appendage 
ere  taken  into  account.  Two  approaches  are  formulated.  The  rigorous 
method  leads  to  a  set  of  two  simultaneous  Integral  equations  which  are 
reduced  to  a  single  Integral  equation  Involving  elliptic  Integrals. 

The  simplified  method  employs  the  technique  of  substitution  vortices 
which  yields  explicit  analytic  expressions  for  the  quasi-state,  apparent- 
mass  and  wake  forces  for  both  the  propeller  and  appendage.  Graphs  and 
tables  give  the  magnitude  and  variation  of  the  propeller  and  appendage 
forces  as  functions  of  tip  clearance,  distance  and  relative  size  of 
appendage  and  propeller- 
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II  INTRODUCTION 

The  purpose  of  this  study  is  to  evaluate  the  effects  due 
to  a  propeller  rotating  in  the  vicinity  of  an  appendage  both  subject 
to  a  uniform  stream  parallel  to  the  appendage,  as  shown  in  Fig.  1. 


. .  ^ 

Figure  1. 

Propeller  Appendage  System 


Specifically  it  is  desired  to  analyze  the  nature  and  magnitude  of  the 
forces  acting  on  both  the  propeller  and  appendage.  These  forces  due  to 
the  rotation  of  the  propeller  will  be  time  dependent  functions  and 
as  such  will  induce  vibrations  either  by  direct  excitation  of  the 
appendage  or  by  the  transmission  of  forces  via  the  propeller 
shaft  The  knowledge  of  the  relation  between  these  forces  and 
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such  design  parameters  at  tip  clearnace,  relative  size  of  appendage 
and  propeller  and  the  location  and  number  of  blades  may  provide  a 
means  for  new  and  more  effective  design  procedures. 

Of  the  early  work  dealing  with  forces  Induced  by  rotating 

1  2 

propellers  are  the  papers  by  Lewis  ’  who  by  use  of  models  and 
electrical  analogy  charted  the  pressure  field  Induced  on  a  wall 
parallel  to  the  propeller  axis.  However  the  results  while  setting 
down  the  nature  of  the  vibratory  forces  do  not  correlate  these  with 
the  various  design  parameters  Inherent  In  propeller- appendage  assem- 

3 

biles.  The  most  recent  work  on  the  subject  Is  that  of  Breslln  ,  who 
solved  the  problem  analytically  by  representing  the  propeller  blade 
by  a  single  concentrated  vortex  sweeping  by  a  moving  plate.  The 
present  work  follows  Breslln* s  approach  but  differs  from  it  and 
extends  it  in  several  important  respects.  Instead  of  a  concentrated 
single  vortex  the  blade  is  here  represented  by  a  chord-wise  vorticity 
distribution  which  in  addition  to  yielding  a  more  correct  value  for 
the  forces  on  the  appendage,  enables  us  also  to  calculate  the  forces 
on  the  propeller  which  was  Impossible  to  do  in  Breslln' s  case.  The 
second  important  difference  is  that  we  have  here  included  the  mutual 
Interference  effects  of  propeller  blade  and  appendage;  for  as  the 
propeller  moves  through  the  velocity  field  of  the  appendage  these 
velocities  modify  and  are  in  turn  themselves  modified  by  the  velocity 
potential  of  the  travelling  propeller. 

Two  approaches  will  be  formulated,  one  rigorous  and  a 
simplified  method.  The  rigorous  analysis  leads  to  a  set  of  two 
simultaneous  integral  equations  which  by  proper  inversion  and 
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Integration  are  eventually  reduced  to  a  single  non-homogeneous  Integral 
equation.  This  equation,  given  In  terms  of  elliptic  Integrals  with 
complex  moduli  and  parameters,  can  be  solved  only  by  numerical  methods 
and  represents  a  formidable  computational  task.  The  simplified  method 
is  the  Method  of  Substitution  Vortices  which  leads  directly  to  explicit 
analytic  expressions  for  the  vortlcltles  and  forces  Involved.  These 
expressions  are  then  evaluated  for  a  number  of  pertinent  parameters 
for  both  the  propeller  and  appendage.  A  discussion  of  these  results 
and  the  main  conclusions  are  given  at  the  end  of  the  report. 

The  analysis  Is  based  on  the  usual  assumptions  of  potential 
theory,  namely  that  the  flow  Is  Incompressible  and  Irrotatlonal  and 
on  the  assumption  that  the  velocity  V  »  U  so  that  the  flow  Is  nearly 
parallel  for  both  the  appendage  and  the  propeller.  The  propeller  is 
treated  as  one  of  a  very  high  aspect  ratio  so  that  its  three-dimensional 
configuration  can  be  replaced  by  a  representative  blade  cross  section 
of  finite  chord.  The  appendage  too  though  Infinitely  long  has  a 
finite  width  in  the  plane  of  the  propeller  chord.  The  problem  thus 
Is  essentially  reduced  to  an  analysis  of  the  unsteady  flow  field 
round  two  flat  plates. 
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The  basic  relations  for  circulation  and  lift  on  a 
flat  plate  are  well  known  from  elementary  aerodynamics.  The 
infinite- span,  inf initely-thin  plate  of  chord  2a  is  represented 
by  an  unknown  chordwise  vorticity  distribution  (^(x)  which 
gives  the  velocity  field  around  the  plate  By  imposing  the 
conditions  of  zero  normal  velocity  along  the  plate  and  the 
Kutta  condition  of  zero  velocity  at  the  trailing  edge  the  velocity 
u(x)  along  the  plate  is  obtained.  Used  in  Bernoulli's  equation 
these  velocities  will  provide  the  pressure  distribution  on  the 
plate  and  by  Integration  these  pressures  will  yield  the  total 
lift. 


F 


Figure  2 

Vorticity,  Circulation  and  Lift  on  a  Flat  Plate 


When  the  plate  Is  subject  to  an  unsteady  free^ stream  velocity, 
a  varying  angle  of  attack,  or  Is  otherwise  affected  by  unsteady 
phenomena,  the  vortlclty  the  i  -  uicn  . ?.nd  tiie 
resultant  forces  become  time  dependent  functions..  Under  these 

4 

conditions  vonKarman  and  Sears  have  shown  that  the  total  force 
on  the  plate  Is  made  up  of  three  parts.  These  are: 

1.  Fq  -  Quasi-steady  Force  This  is  the  force  that 
would  have  resulted  from  a  steady  state  case  in  which  all  para¬ 
meters  are  equal  to  the  Instantaneous  values  of  the  unsteady 
motion.  The  expression  for  this  force  is 

-  pu  r„  (3-1) 

2.  "  Apparent-Mass  Force.  This  is  the  inertia 
force  on  the  plate  due  to  the  acceleration  of  the  mass  of  fluid 
around  the  plate.  It  is  given  by 

-  -  P  ^  /a  X7^(x,t)dx  (3-2) 

3.  ^2  "  Force,  Ihie  to  the  variation  of  7^(x,t) 

with  time  the  foil  will  shed  vortices  behind  it  giving  rise  to  a 
vortlclty  distribution  7(C,t)  the  wake.  With  respect  to  the 
coordinate  axes  of  Figure  3  this  will  be  a  function  of 

both  t  and  ^  although  in  space  it  would  not  be  a  function  of 
time.  This  new  vortlclty  distribution  will  give  rise  to  a 
normal  velocity  component  on  the  foil  which  must  be  supressed  by 
an  additional  vortlclty  72(x,t)  on  the  foil. 


Figure  3. 

Wake  Vorticity  Due  to  T"  nnstu;'.’'  Ci.rovjlation 


It  la  shown  in  Section  VOL  that  the  circulation  in¬ 
duced  on  a  plate  2a  by  a  single  vortex  at  z  is  given  by 


For  the  case  of  a  single  vortex  located  at  y~0  this  becomes 


r'(x,o) 


1 /x-a 
V  x+a 


1) 


For  a  vorticity  distribution  7(^,t)  along  the  wake  we  thus  have 


where  ^  R  "  And  r2  Is  the  additional  circulation  on 

the  plate  due  to  the  wake  vorticity. 

We  now  have  three  circulations  to  consider;  two  about 
the  foil  given  by  and  and  one  in  the  wake,  P,  caused  by 
the  unsteady  motion  of  the  foil.  By  Kelvin' s  theorem  the  sum 
of  these  three  circulations  must  at  any  time  be  zero.  We  thus 
have 


r  +  p,  +  I'  =  0 

o  2 


or 


1+R  _  /T  1  _  1+R  __ 

r„  +  a  /  t'(C.t)  (  /^-Ddc  +  a  /  /(C,t)dC  -  0 


C+1 


which  yields 


f  ^  dc  =  -  -2—  (3-4) 

1  V  ^+1 

The  above  is  an  integral  equation  in  y(^,  t)  which  must  be 
solved  for  the  particular  function  I^(t)  given  by  (2-11). 

The  force  on  the  foil  due  to  the  additional  circula¬ 
tion  r„  can  be  shown  to  be  given  by 


1^-R 

F  ^  fUa  J  /(^,t) 

1 

The  total  force  on  the  foil  is  then 


AL 


/n 


(3-5) 


F  =  F  +  F,  -»■  F, 
o  1  2 


(3-6) 
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IV.  FORCES  ON  A  SINGLE  FOIL 

In  this  section  the  vortlclty  distribution,  velocity 
potential  and  the  forces  arising  from  the  iiiot*. on  oi  c 
single  foil  in  a  perfect  fluid  will  be  set  down  to  serve  as  a 
starting  point  for  the  treatment  of  the  more  complex  propeller- 
aV’pendage  problem  The  results  derived  here  will  also  be  needed 
later  on  in  the  text  The  elementary  model  considered  is  that 
■- 1  a  simple  flat  place  moving  with  a  velocity  U  at  some  attack 
angle  a  aa  shown  in  Figure  4.  In  this  sketch  we  have  the 
coordinate  axes  fixed  with  respect  to  the  plate  and  we  let  the 
jrream  flow  past  the  foil  with  a  velocity  Ue^*^, 

I’ 


Cl  I  OL 

i 


Figure  4 . 

Fiat  Plate  in  Uniform  Stream 
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The  complex  potential  for  a  simple  concentrated 
vortex  of  strength  located  as  z  ■  Is 

^0 

and  for  the  uniform  stream  U 
w  ■  -  Ue^®z 

From  the  above,  representing  the  flat  plate  by  a 
distributed  vortlclty  y^Cx),  we  have 

w  -  -  Ue^“z  +  ^  yQ(x')fn(z-x’)dx'  (4-1) 


From  (4-1)  the  total  compiej;  'clocicy  ii  th«(. 


dw  _ 

3?  “ 


Ue^®  +  f 
Zv  - 


a 


a 


dw 

^  -  u  +  Iv  we  obtain 


Since  ^ 
and  normal  to  the  foil 


7^(x')dx' 


for  the  velocity  components  parallel 


u 


Ucos  a 


®  7„(x')dx' 

r  o 

(,:x’)2+y2 


(4-2) 


1  ®  7o(*’) (x-x') 

V  -  Usln  a  +  ^  /  - - -T-^—  (^-3) 

(x-x')  +y 

The  condition  of  zero  normal  flow  along  the  plate 
requires  that 

v"0  at  y-0  -a^x^a 


and  thus 
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1  7o(x')dx'  _ 

_  j  I  " 

TT  -a  x*x' 


2  Usln  a 


(4-4) 


This  represents  an  integral  equation  in  y^i^)  and  its 
solution  can  be  obtained^ from  the  following  general  inversion 
formula.  If 


i  /  -  v(z)  (4-6) 

TT  S-2 

®1 

then  together  with  the  Kutta  condition  v(a2)  =  0,  we  have 


For  the  case  of  equ  (4-4)  we  thus  have 


/  *  oj-ii  ifa-v\  r  I  / a~^'X  dx 

7^(x)  » - ;; -  J  y-T^  p-— 


2U  sin  a  Ja-x 


Since 


(4-7'< 


^  ^  -d2s:  .  , 

^  Va-x'  x'-x 

a 


(4-8) 


we  have 


7o(x) 


2U  sin 


“  \/l^x 


V  4  ' 


For  non-steady  flow  the  right-hand  side  ot  eq 
is  given  by  an  arbitrary  function  of  space  and  time 
Hence  the  vorticlty  distribution  along  the  foil  will  also 
a  function  of  time  and  the  Inversion  integral  then  beeoa 


The  total  circulation  around  the  foil  is  then 

a 

r^Ct)  -  /  7Q(x,t)dx 


As  seen  from  the  above  simple  example,  the  . i 
in  arriving  at  an  expression  for  the  aerodjmamic  foiv 
represent  the  foils  by  some  vorticlty  distribution  y(/ 
into  account  all  neighboring  velocity  potentials,  Ih^ 
applying  the  proper  boundary  conditions  the  unknown  v. 
distributions  can  in  principle  be  obtained  from  the  at 
relations.  These  procedures  will  now  be  applied  to  th 
on  hand. 


V.  RIGOROUS  ANALYSIS  OF  PROPELLER- APPENDAGE  SYSTEM 


We  are  concerned  here  with  the  unsteady  potential 
generated  by  a  propeller  rotating  in  the  vicinity  of  a  moving 
plane  or  appendage,  as  shown  schematically  in  Fig.  1.  The 
horizontal  distance  of  the  blade  from  the  appendage  is  a 
fixed  quantity  which  can  be  either  negative  or  positive  while 
the  vertical  distance  y^  is  a  variable.  Thus  the  velocities 
and  forces  generated  on  both  the  appendage  and  propeller  will 
be  time  dependent  functions „  A  crucial  element  of  the  problem 
is  that  there  is  mutual  interference  between  propeller  and 
appendage.  In  addition,  due  to  the  unsteady  circulation, 
wakes  will  be  generated  behind  both  propeller  and  appendage. 

Mathematically  the  system  can  be  represented  by  two 
flat  plates  at  right  angles  to  and  at  a  continuously  vatfing 
distance  from  each  other,  as  shown  in  Fig.  5.  We  shall  use 


Figure  5. 

Schematic  Diagram  of  Propeller  Appendage  System 
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the  following  notation  for  the  various  vorticities  involved: 
7p^(y)  -  vorticity  distribution  on  propeller 
at  infinity 

7pi(y)  -  vorticity  distribution  on  propeller  due 
to  appendage  interference 
7op(y)  ”  quasi-steady  vorticity  on  propeller 
7'  (x)  -  vorticity  distribution  on  appendage  due 

Ocl 

to  7pi(y) 

“  vorticity  distribution  on  appendage  due 

Oa 

to  7p„(y) 

7  (x)  -  quasi- steady  vorticity  on  appendage. 

0& 


With  the  coordinate  system  fixed  in  the  center  of  the  appendage  and 
the  main  flow  parallel  to  x,  i.e.  a  ■>  0,  we  have  from  equ.  (4-1) 

<  1 

7iF  (  7op(y')in(z-x^j-iy')dy' 
i  ^ 

(5-1) 


The  complex  velocity  potential  is  by  differentiating  (5-1). 


dw 

o 

ST 


7oa<*'>dx’ 


7op(y')dy' 

z--tx;+iy'') 


or  after  rationalizing  the  denominators 
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.  [(x-x  )-i(y-y') ] 

+  7^_(y')  - 9 - 9—  dy’ 


The  component  velocities  are  thus 


u  -  U  -  11- 
Zir 


TT 


^  (x-xO^'V  ^ 


.a. 


(x-xj  +(y-y') 


v-f 

Ztt 


,  (x-x')^+y^ 


)®2  (x-XQ)dy' 


(5-2b) 


The  boundary  conditions  to  be  satisfied  are 

1)  On  appendage:  v«0  on  y“0  ^  x  ^  t 

2)  On  propeller:  u"0  on  x-x^  yo 

Thus 


224^0^ 

x’  -X 


2  7op(y')(x-x^)dy' 


(5-3a) 


U 


7^ 


7o-(x')dx' 

y  /  — 72^  + 

-/  (Xo’*'>  +y 


a 


a. 


2  7op(y')dy' 

y-y' 


(5- 3b) 


Since  7op^y'^  "  ^poo^y*^  "^pl^y*^  rewrite  the  last  Integrals 

of  the  above  two  equations  as  follows: 
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(x-Xo) 


A  7op(y’)dy' 

Ja^  (x-Xp)^+(y')^ 


(x-Xq) 


fl  y  Ay') dy* 

3^  i^-^Hy')^ 


+ 


7p«>(y')dy* 
(x-x^)  +(y') 


1  A  7pp(y')dy’  _  1  ^2  7pi(y')dy'  ^  /?2  7p„(y’)dy’ 

-Vf—  -  ^  J  ^  /  y^y^ 

u  (/  a^  0  a^ 

The  last  Integral  is  of  the  form  given  in  (4-4)  where  Usin  a  represents 
the  cross  flow,  in  our  case  merely  U.  Thus 


1 

I? 


(y')dy' 

yTyH 


1 


u 


(5-4) 


Hence  in  Equ.  (5- 3b)  the  term  U  will  be  cancelled  by  the  7p5g(y') 
portion  of  the  second  right  hand  integral.  Furthermore,  from 
Equ.  (4-9)  the  value  of  7p„(y*)  satisfying  (5-4)  is 


7  (y)  -  2  U 


(5-5) 


Thus  after  introducing  the  value  of  7p„(y')  into  (5-3a)  a  portion  of 
that  integral  will  become  after  integration  a  known  function  of  x 
alone.  We  can  now  rewrite  (5-3)  as  follows 
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S„(x')dx' 
.  x'  -X 


®2 

(X-Xg)  / 


7p£(y')dy' 

(x-x^)^+(y')^ 


+  f(x) 


(5-6a) 


/  'yoa(x')dx' 

y  J - 5 — y 

-I  (Xjj-x')^+y^ 


2  7, 


(y’)dy' 

y'-y 


(5- 6b) 


Employing  the  inversion  formula  (^-7)  for  the  above  two  Incegral 
equations  with  the  Kutta  condition  applied  at  |  for  the  appendage 
and  at  ^2  the  propeller,  we  have 


'”■■7  l/v?  i 


Vy- 


a,-y  y’-y 


J  >^(x’)dx' 

/  -  - - ^  dy' 

-i  (x„-x)+(y') 


(5-7) 


y'i 


6a' 


(x)  -  - 


rtz£ 

i+x 


Es  (•-*o> 


\/^ 


s  s-x 


®2  [7pi(y')+7p„(y')]dy' 

J  - - - C - - - 


ds  (5-8) 


Breaking  up  the  appendage  vorticity 7^(x)  into  ay^^  (x)  corres¬ 
ponding  to  7pj^(y),'  and  a-/^^  (x)  corresponding  to  7p„(y)  we  can 
rewrite  (5-8)  as  follows 
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i 

s 

-i 


8-X 

( 

8-X 


t 

s 

-t 


8-X 

_ o 

8-X 


/  - - - 7 

®1 


d8 


®2 

/ 

®1 


^a.(y*)^y 

(s-Xq)  +(y') 


'  d8 


(5-9) 


(5-10) 


Integrating  first  with  re8pect  to  s  we  are  confronted  with  the 
following  Integral 


I 


t 

S 


da 

- ? - 7 

(8-x^)^+y'^ 


(5-11) 


This  Integral  which  reappears  several  times  later  on  Is  solved  In 

general  terms  In  the  Appendix.  The  solution  for  the  case  of  com- 

2  2 

plex  conjugate  roots  of  the  quadratic  (s-x^)  +y'  Is  given  by 
Equ.  (A-9) .  For  our  case 


ki  -  - 


1,2 


*0  -  iy' 


Thus 


I 


V  Re 


1 

(Xo-x)+iy' 


(x^+f)+iy’ 

(Xp-l)+ly' 


(5-12) 


and 


7Afl(*) 

oa 


i  \/ilx 

TT  V  l+X 


a. 


Re  / 


a. 


(Xjj-x)+ly’ 


(x  +i)+ly' 

T^r-ipiy'  7pi(y')dy'  (5-13) 
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(5-14) 


Evaluation  of  71 • (x) 

- - -  oa - 

Introducing  the  value  of  7p^(y)  as  given  in  Equ.  (5-5)  into  (5-14) 
we  have 


where  =  i(x^+0  b2  =  i(Xp-^)  b  =  i(x^-x) 


The  integration  of 


is  performed  in  the  Appendix  and  its  solution  in  terms  of  elliptic 
integrals  is  given  by  equ.  (A- 15)  as 


s  t 

where  K(k)  is  the  complete  elliptic  Integral  of  the  1  kind, 
|T(e^,k)  the  complete  elliptic  Integral  of  the  second  kind  and 
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V  ±m  '  /  •  ^ 


^11^12 


.  yo-i(V« 

a,,  - - 

7 


‘11 


a 


‘13 


c 

7 


_  yo-i(v^> 

12  - 5 - 

7 

(1-ailv) 


1-va, 


.  (IJZli.) 


ai3-v 


The  explicit  expression  for  then  becomes 


=  ity  >/iEE 

Im  l'*’F 

,1  f, 

TT  Y  £+x 

"*  v(l-va33) 

V  ^®ll‘F)(a32*v) 

(5* 


+  (l-v)(Uva^3)  /Tcv^.k)  ^  ai3)(l+ai3) 


F.valuation  of  7pj^(y) 

The  expression  for  7pi(y)  is  given  by  Rqu.  (5-7).  Since 
/ygCx)  is  a  function  of  x  only  we  can  rewrite  it  as  follows; 

vo-)  ■  -  7  {W,  {[  i 


■'  '  dx' 

>"  (X  -X')'  +  (y.)^ 


17) 
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y  •  -y 

Writing  5  the  integral  in  y'  reads 

C 

7 


with  2  "  “  I 


We  see  that  this  integral  is  of  the  same  form  as  (5**  11)  and  its  solu¬ 
tion  for  complex  conjugate  roots  is  according  to  (A-9) 


I 


Evaluation  of  7i„(x) 

*oa 

Introducing  now  Equ.  (5-18)  into  (5-14)  we  have  for 
the  following 
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Since  Re  Re  Z2  ■  -j  (Re  ^-^^2  +  Re  have  after  multiplying 

top  and  bottom  by  -i  and  writing  Im  z"Re(-i2). 


In  abbreviated  form  the  above  reads 
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It  It 
^al  FTg  Ib2 

where  - (’:^+i)  b2"i(x^-|) 

b  -i(x^j-x)  g  -i(x^-x') 

Integrals  1^2  ^^2  same  form  as  the  Integrals 

of  F.qu.  (5-15)  and  their  solution  is  given  by  Equ.  (A- 15)  with 

'"’U  “  I  ^  replacing  in  1^2 

®15  “  c  ^  replacing  a^^^  in  I^j2 

We  thus  have 
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T  '  I  ^  (l-a,,v)K(lc)  + 

^al  (l-vaT^T 


(l-v)(l-va. 


^a2  "  V  (aii-v)Ui2“^T 


;  (l-a^^v)K(k)  + 


L 


(1-v)  (1-va^^)  II  (v  ,10  +  * 


T  .  2(lj30 

^b2  (l-va^jl 


l»v 


(aii-v;(a^2“^ 


(l-a^^v)K(k)  + 


(1-v)  (l-va^j)  IF  (-/^k)  + 


v(l”v) 


ai5-v 


n  (8^.io| 


the  expressions  for  and  1^  and  then  In  Equ.  (5-19),  we 


Using  these  in 
have  for  the  following: 


“  ■T  \/M  \/laji-viia^2^ 


i 

-i 


a. 


^  iFo^.lO-kjlTcs^plOldx' 


a2-lix^-x^y 


\  IksK(k)  tk2"f[  0^,k)-k3ll (7^,10  ]dx' 

\|a2+l(xQ-x’T  x’-(2x^-x7  '  1  '  Z 


(5-20) 
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Functionally  Equ.  (5-21)  can  be  written  as 


yf  /X'i  -  -ir  V  /i.Ti£ 


[7ia<*’)+7’t;a(x’)l  f(x',x)dx' 


which  is  a  non-horaogeneous  integral  equation  in  7Qa(x)  with  7Qa(x) 
given  by  Equ.  (5-17). 
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V I  SIMPLIFIED  ANALYSIS  -  METHOD  OF  SUBSTITUTION  VORTICES 

It  is  apparent  from  the  above  that  the  rigorous  approach 
yields  analytically  unmanageable  expressions.  If  at  all  feasible;, 
the  above  integral  equations  would  have  to  be  solved  numerically 
which  because  of  the  finite  difference  expressions  and  iteration 
procedures  in'olved  would  appreciably  affect  the  accuracy  of  the 
results  We  shaLs  therefore  use  a  different  approach  to  the 
problem  which  ol  though  invcl'^ing  certain  initial  simplifications 
permits  us  tc  obtain  explicit  analytical  solutions  with  an 
accuracy  perhaps  higher  than  would  have  resulted  from  a  numerical 
evaluation  of  the  exacu^  equations  While  bcth  the  appendage 
and  the  prcpeller  are  treated  as  finite  length  plates  with  dis¬ 
tributed  vcrricities  the  effect  that  each  of  these  plates  has 
uC'  its  neighbor  is  assumed  tc  be  that  of  a  concentrated  vortex 
whos®  strength  equals  the  total  distributed  strength  of  the 
respective  plate  The  -pethod  of  evaluating  the  unknown  strengths 
of  the  propeller  a-'d  appendage  will  become  clear  in  the  course 


of  the  following  analysis 

AS  explained  previously 


the  propeller  experiences  at 
infinity  due  to  the  cross 
flow  U  a  vorticlty  y) 

given,  hy  equation  (4  9)  As 
the  propeller  approaches  the 
appendage  it  begin'^  tc  feel 


oa 


Fig  '  Representation  oi  ,  (^x)  by 

v-c*. 

a  concentrated  vortex  F, 

'  i-t 
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the  effect  of  7-.(x) ,  chus  giving  rise  to  an  additional  vorticity 

Oct 

7pi(y) .  Referring  to  Fig.  6  we  shall  now  concentrate  the 
distributed  vorticity  7..(x)  of  the  appendage  at  some  location  z 

Ocl  o 

and  call  it  r  .  Since  7  cancels  the  effect  of  the  uniform 
oa  poo 

stream  U  we  have  for  the  complex  potential  due  to  a  concentrated 
vortex  and  a  distributed  vorticity  7p^  the  following: 


ir. 


0^2 


W  -  In(z-z^)  +  ^  'Vpi(y')  in(z-Xp-iy')dy’  (6-1) 


a. 


From  this 


dw 


27r(z-z^) 


1  y.Ay') 


Applying  the  condition 


u  •  0  at  X  ■  x^  ^  y  ^  ®2 


T-y 


r 

^  oa  Z  \  2  ,  77 

(x-x^)  +(y-y^) 


This  Integral  equation  according  to  (4-7)  yields 


s  (s-y^)^+(Xj,-x^)^ 


ds 

s-y 
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The  circulation  is  given  by 


a. 


rq,  -  ( V  "'pi’  dy  -  2U  / 


y 

y  +  f  -  Vo 


dy 


r  ^2 

’ex; 

J 

TT 


yo+ 1 


■yo  +  7 


y-Yo  +  7  y  Yo  +  7  <®‘ya^ 


ds 

7  s-y 


y_y 

Using  the  substitution  s  *  - ^  we  have  for  the  first  right-hand 

•)» 

/ 

integral 


2U  I  f 


Y+f  dC  -  ^UC 


(6-2) 


For  the  second  integral,  integrating  first  with  respect  to  y 


7 

this  being  a  Cauchy  princiiJdL  value  integral.  Now  since 


„  1 

- ^ - .  Re - 

Is-Ya)  +(>o-Xa/  [y^+Kx^-Xa)  ] 


and  writing 

C  - 


s-y  (y  -y  )+i(x  -x  ) 

_ ;^_u  ,  ^  ■'a _ ^  o _ a;_ 

c  ^  c 


we  have 


r ,  1 

Re  S 
-1 


i+c  dr 


1 1  a  ■! 

■  r 


(6-3) 
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The  solution  to  the  above  Integral  is  given  by  (A- 7)  and  thus 

^pi  "  ■  ^oa  "  x/f^T  ^ 

The  total  circulation  is  thus 


op  p«o  pi 


ttUc  +  r  Re  [ 
oa 


<vya" 


-  1  ] 


(6-4) 


To  find  7Qa(*^  treat  7Qp(y)  as  a  concentrated  vortex  located  at 
some  appropriate  coordinate  and  write  the  complex  potential  as 
due  to  r^p  and  a  distributed  vorticity  7Qa(x)  along  the  appendage. 
We  thus  have 


„  .  in(z-x  -ly  )  +  ^ 


7Qg(x')  in(z-x')dx' 


Figure  7. 

Representation  of  7Qp(y)  by  a  concentrated  vortex  r^p 
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dw  -  1  ,  i  / 

a?  -j^r  J-’Jp'iy;;  — 


The  condition  of  no  normal  velocity  along  the  appendage  requires 


and  thus 


V  -  0  on  y  -  0  x 


I  7  (x')dx'  r  (x-x  ) 

r  ut.i  p 

j  -  I"  *  .■"  *  ■  "TT 

(k-x7^v 


Again  using  F,qu.(4-7)  for  this  integral  equation  we  obtain 


7  U/ 
oa 


,/IS  /  ,/Ei 

7  v*+x  Z,  Vi-S  (3  ) 


)  +y 

p 


(6-5) 


Since 


S  (x)dx 


oi  oa 


we  can  write 


r  * 

r  -  -  -5?  / 

oa  ^  , 

TT  -I 


flEE  ^  s  Jm 

f+X  S-X  _  .  V  i”8 


(®-V 


' 

TT 


/  ,/E| 


By  partial  fractions  the  above  integral  can  be  broken  up  into  a 
sum  of  two  integrals  of  the  form  of  I,,  with  their  solution 
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where 

B  - 

The  method  of  determining  x  ,  y  ,  x  and  y  the  coordinates  of  the 
substitution  vortices^ is  given  in  detail  in  the  Appendix.  The 
exact  formulae  for  the  four  coordinates  are  given  by  equations 
(B.15),  (B.18),  (B.20)  and  (6-8).  It  can  be  verified  either 
analytically  or  by  sample  calculations  that  x  ■  x  and  y.  “  0  ,  i.e. 

^  w  O 

that  the  substitution  vortices  lie  on  the  plates,  the  only  variation 
occurring  along  and  not  away  from  the  appendage  or  propeller.  We 
can  thus  rewrite  A  and  B  as 


<yo-  f)^^<^o-^a> 
(yo+  7)+i(Xo'*a> 


(6-9) 
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(7-4) 


As  mentioned  in  Part  III  the  determination  of  the  apparent- 
mass  forces  j.nvolves  the  evaluation  of  the  moment- of-momentum 
integrals  namely 
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Using  for  7^-(x)  the  expression  (6-5)  we  have 

Oct 


\l  £-s 


S  0 


(s-Xp)  +yp2 


The  s  integral  above  has  been  solved  previously,  equ.  (A-1),  thus 


M  »  _2£  Re 
a  IT 


[X  +i)+iy. 


J., 


x„-x)+ly. 


Let  us  consider  the  above  integral  with  z  “  x  +iy_ 

H  r  r 


■  />  >/ 


^  dx 
i+x 


(7-5) 


Breaking  it  up  by  means  of  partial  fractions 


■  /  \l^  ==1 


ft  , _ 

I  .  /i-x  dx 

1 , 


and  using  for  the  first  integral  the  substitution  x  =  cos  0  and  for 


the  second  the  solution  given  by  equ.  (A- 7)  we  have 


I  =  -  XtT  -f-  ZpTT  [  1- 


Tluis 


=  r  Re  [ 
n  op 
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(7-6) 


(7-7) 


For  Mp  using  the  expression  for  7^p  we  have 


ar 


c 


M  =  2U  /  (y-y^) 

*  fl 


r 

oa  /■  /  X 

■  — /  <y-yo> 


y  +  f  -y^ 


dy 


y^^  2-y 


I 


s-y 


a 


y.y^+l  aj  (s-y^)  +(x^-xp' 


ds 


Writing 


y-y. 


c 

2 


=  ^  we  have 


2  1 

M  =  2U(|)  /  C 


-1 


1+C 


2 


£(s)ds 


-  C 


=  2U(f)"  (-  |)  - 

TT 


r^2 


a. 


7r[(s-y^)  -  I  )  f(s)ds 


,wCv2 

=  -  ttUC^)  +  ~ 


(^2 


(s-y„)^  Re 


o  a 


Using  the  substitution 


c 

2 


^  we  have  for  the  integral 


where 


1=1  /,  ^  =  f  ^t-P'  -^v/^  1 


-1 


(7-8) 


P'  = 


(y  -y  )+i(x  -X  ) 

c 

2 


Thus 


oa 


o  'a'  '  o  a 


o  '  a  o  a 


=  -  TrtJ(§)‘+ r^^Re  Uy«-yJ+i-(x„-xJ- j[(y„-yj+i(x„-xjl‘’-(|)‘  ) 


(7-9) 
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The  apparent  mass  forces  are  given  by 


“P 


A. 

ht 


P 


at 


Since 


at 


V 


ay 

o 


Since  x  =  x  and  y  =  0  we  can  rewrite  the 
p  o  a 

normalized  expressions  for  and  Fj^p  as  follows: 


la  d  - 

~  '  ST  '^op 
o 


r_  Re  [  y(x^+iyp)^-  1 


^  =  ST  ^  } 

Xq  °p 


(7-10) 

(7-11) 


In  differentiating  the  above  expressions  it  should  be 

kept  in  mind  that  these  in  addition  to  y  contain  y  and  , 

o  P  ^ 

themselves  functions  of  y^  as  given  by  equs.  (B.18)  to  (B.20) 
in  Appendix  B. 
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c)  The  Wake  Forces 

As  mentioned  in  part  III  the  unsteady  circulation  on 
both  the  appendage  and  propeller  will  cause  the  shedding  of  wake 
vortices  and  these  in  turn  will  Induce  an  additional  circulation 
round  the  two  plates.  For  the  appendage  the  circulation  caused 
by  a  single  concentrated  vortex  located  at  some  distance  x  along 
the  wake  is  by  setting  c  ■  y^  -  0,  *  x  in  equ.  (6-7)  given  by 

-1)  (7-12) 


Thus  according  to  cqus.  (3-4)  and  (3-1))  the  integral 
equation  for  the  wake  vorticity  y(0  the  force  induced  by  it 
on  the  appendage  are 


(7-13) 


l+R 

F,  -  pU«  / 

^  1 


(7-14) 


Introduce  np,  the  variable  s,  measured  from  the  initial 
position  of  the  appendage, s  «  (l+R)- C, we  have  for  equ  (7-13)  by 
writing  p,(s)  -  7(l+R-s) 


Coordinates  of  Wake  Vorticity 
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R 

/  M-(s) 


0 


(7-15) 


This  Integral  equation  Is  difficult  to  solve,  particularly  In  view 
of  the  complexity  of  the  expression  for  We  shall  obtain  a 

solution  to  (7-15)  for  a  unit  Impulse  l.e.  for  =  1  and  then  sum 
the  result  by  means  of  the  Duhamel  Integral.  For  a  unit  Impulse 
equ.  (7-15)  reads 


I 


(7-16) 


and  thus  the  unit  force  produced  by  this  unit  impulse  Is 


pU)  . 


1 


/ 


(7-17) 


r-1 


where  is  the  solution  to  (7-16). 

According  to  Garrick^  the  value  of  the  Integral  in 
equ.  (7-17)  is  within  an  accuracy  of  27,  given  by 


4+R 


and  thus 


a  -  (  unit  circulation)  (7-18) 

The  above  Is  the  force  due  to  a  unit  Impulse  of  magnitude 

r  aO  at  t  <  0  and  F  =1  for  t  >  0.  Using  the  Duhamel  method  of 
oa  oa  ^ 
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superposition  for  any  arbitrary  input  we  have 
AFj  -  F^^^(e-T)Ar„^<T)  +  r„^(0)F^^\t) . 


dr 


Od. 

Since  at  t«0  T  mQ  and  aT  (t)  =  —tz~  dr 
oa  oa  OT 


2a 


2pU 


^  dT 


Integrating  by  parts 


dr 


4  +  j(t-T) 


(7-19) 


2a 


Writing 


2pU 


^2a  " 


r  .1 

t 

.t 

r 

_  4  .  ^(t-T) 

0  > 

oa 

o 

-  Ut 

"  =T 

[‘2  oa 

(4+R-t) ^  ^ 

d(^T) 


4  +  j(t-T)  ^ 


or 


rF 


2a 


1 

-  ( — 2®)  +  2 
2  V  T  /  + 


aR  (T)dT 

oa  _ 

.  o  (4+R-t) 


(7-20) 


Eq.  (7-20)  tells  us  that  the  wake  force  is  always  equal 
to  -  j  or  minus  one-half  of  the  quasi  steady  force  plus  a 
correction  C  given  by 


C  =  2 

In  the  above 


oa 

o  (4+R-t) 


(7-21) 


f  oa‘^>  “  ye(A-{)ReiB-l) 

where  v  -  a“^^  with  t  replacing  R  for  purposes  of  integration. 

»  r\ 
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The  constant  a  Is  the  Initial  position  of  the  propeller,  later 
on  to  be  pushed  to  Infinity.  In  this  expression  t*0  implies 
that  the  propeller  is  at  •  a  as  shown  in  Figure  9.  In  order 
to  have  the  propeller  at  y*0  at  t"0  we  shall  employ  a  new  time 
variable  to  be  related  to  t  by  the  expression 

R  -  \  I 

with 


Using  chis  in  (7-21)  we  have 


Figure  9 

Time  Coordinates  for  Propeller  Position 
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Direct:  Integration  of  (7-22)  is  difficult  since 
in  the  numerator  is  a  complicated  function.  However  we  can 
represent  by  a  simpler  function  which  by  a  suitable  choice 
of  constants  can  be  made  to  give  a  fair  approximation  to  the 
original  expression. 

r  in  general  has  the  shape  of  Figure  10  and  we  can 
Od 

write 


r 

oa 


q  r 

^a  am 


(y  -y  )  +q 

o  om  ‘a 


(7-23) 


where  q  is  a  constant  so  chosen  as  to  best  approximate  the 

exact  r  .  Since  for  large  y  T  ~ const.  ,  r  F  Re'A-1) 
oa  °  'o  op  oa  op  ' 

which  when  expanded  and  restricted  to  terms  no  higher  than 
yields  (for  -  f )  ° 


Figure  10. 
Shape  of  r  (y  ) 
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Thus  for  large  we  have  q.r _  ■  -  -i  and 

”  o  ^a  am  o  L 


am 


Introducing  (7-23)  and  (7-22)  and  remembering  that 
rr 


R-  +  ^ 


2  T  / 

^a  am 


zr~Y 


(4+Rq+  f  -  t) 


To  eliminate  a  from  under  the  Integral  sign  write 
■'  om 

and  we  have  after  letting  a  <» 


C  -  2rq^r  / 
^a  am 


dt 


<Vyom> 


(7- 


where  k  -  4r  -  (y^- y^^) 


Breaking  up  (7-24)  into  partial  fractions  we  have 


C  - 


2q  rr 
^a  am 

0?^ 


-k  / 


2rdr  ^  ,,z 


p  r+q. 


P  r+q. 


00  00 

+  (R^+q.)  /  — n  +  2k  / 

®  P  (C+k)^  P 


where  P  ■  (yo*yom^  These  Integrations  are  easily  performed 
yielding 
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[(4r-P)%^]^ 


( 


(4r-P)  +q 

(4r-P)  in - 5—2 

16r 


(4r-P)-^.q 


a 


(I  .  tan-l  + 


(4r-P)^+q,  ) 


(7-25) 


Thus  the  force  on  the  appendage  due  to  the  wake  Is 
— “  7  ^ ^ 


(7-26) 


For  the  propeller  the  effect  of  an  arbitrary  concentrated 
vortex  located  at  some  (x,y)  Is 


Re  ( 

op  O  ' 


(yo-y-^+t(x<,-^  _ 
(y„-y''^+i(x„-x) 


Shifting  the  origin  to  we  obtain 


^nn  -  Re  (  \  /  ^ 

°p  °  V 


C„  Re  (  \/  .1  ) 

c+iq 


t 


4^ 


Figure  11. 


Propeller  Coordinates 
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In  accordance  with  our  original  assvmptlon  of  V  »  U 
the  shed  vortices  In  the  wake  are  assumed  to  lie  essentially 
along  T)  >  0 ;  thus 


L„(C»o)  -  C  (  -  1) 


op 


C-c 


This  Is  an  expression  similar  to  (7-12)  and  we  can  thus  write 
down  the  expression  for  the  wake  force  as 


1+S 


Fap  -  pu  (f)  /  ^ 

K  -1 


(7-27) 


where  S  - 


and  thus 


Vt 
^  • 

The  treatment  of  (7-27)  Is  the  sane  as  that  of  (7-14) 


‘  4  F  +  ZirpUVc  /  -32 - — 

2  op  ^  -'o 

Vt  ti 

where  t'  ■  and  all  quantities  In  have  been  normalized  by 
c/2  Instead  of  by  i,  and 

y  '  »  a'  -  t'  . 

■'o 

Using  S  -t-  a'^a  shift  of  variables  to  a  new  time 

t^  clmilar  to  that  illustrated  In  Figure  9^ we  are  now  to  evaluate 

S  +a‘ 


r  (f) 

-32.  ,  dT' 

V 


'o  (4+S^+a  -t') 


TECHNICAL  RESEARCH  GROUP 


49. 


The  shape  of  r^p 


(v^)  as  shown  later  on  is  of  the  form  given  in 


Thus  r  _  -  +  r  . 

op  poo  pi 


and 


S  +a*_ 

°  r  +  r  . 


'2p  -  -  ?  ^op  " 


r  _ES2 - Ei--_  dT 

J  (4+S+a'-T)^ 


S  +a'  . 

O  p 

i  fop  ^  i  fpco  *  2a/>UVc  -£i - 

V  (4+S  +a*-T' 

o  o 


df 


)‘ 


or 


S  +a* 
o  . 


!2e 

T 


-  I  fopl/f  ^  2 


(4+Sjj+a'-'r') 


dx 


Thus  since  Tp^^  is  the  only 
total  circulation  the  wake  force 
correction  terra  given  by 


varying  compouant  of 
is  equal  to  -  ^,)pi 


the 

plus  a 
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S  +a' 
o 


D  -  2  / 
o 

In  the  above 


(^S^+a 


TIT? 


dx' 


O 


-  .  X  RerA-nRe(B-l) 

l-La-ljAeTB-TT 

where  jfor  large  values  of  Yq  can  be  shown  to  vary  as 

c(Xo- 

— 


To  eliminate  a*  from  under  the  integral  sign  we  write 
a'-T*  =  i 

and  since  - 

o  o 


D  =  -  2  f° 


a'  (4-y^  + 


Letting  a'  — ^oo  we  have 
0=2 


/»  1  .  ( 

[ 


(^-y^  +  0‘ 


-28) 
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VIII  NUMERICAL  EVALUATION  OF  FORCES 

By  sumnlng  the  three  component  forces  we  have  for  the 
resultant  force  on  the  appendage 


or 


rF„  1  a 

_ a  _  1  /_oa\  .  la 

"Tf”  ■  7  ( — W")  + 


+  C 


(8-1) 


This  expression  Is  Implicitly  dependent  on  r  and  Is  then  a  function 
of  four  parameters  y^,  x^,  c  and  r.  The  total  force  on  the  propeller 
is 

F  F  F,  1  F  . 

^  j +  D 

or  since  F„p  -  Fp„  +  -  T  +  F^j 


^-1.0  +  ^!^  +  ^+D  (8-2) 

The  force  on  the  propeller  due  to  neglecting  second  order  effects 
is  Independent  of  r. 

Equations  (8-1)  and  (8-2)  were  evaluated  on  an  IBM  7090 
digital  computer  using  the  expressions  derived  In  Section  VII. 
Expression  D  as  given  by  Eq.  (7-29)  was  the  only  unevaluated  integeral 
and  It  has  been  Integrated  nxnnerlcally  using  Simpson's  rule.  The 
following  discrete  values  of  the  Involved  parameters  were  used 
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y  from  +  5.0  to  -5.0  In  Intervals  of  0.1 

o 

c  0,  .Ij  .2,  .3,  .5 

\  ±l-6»  ±1-4,  +1.2  +  1.05 

r  2,  3,  4. 

The  values  for  both  the  component  and  total  forces  are  available 
In  tabular  form  and  the  essential  results  of  these  calculations  are 
plotted  In  Figures  13  through  51.  In  these  all  forces  are  plotted 
In  terms  of  the  unperturbed  propeller  thrust  T^or  whereas  all 

distances  and  the  propeller  half-chord  c/2  have  been  normalized  by 
the  appendage  half -chord  i. 
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IX.  DISCUSSION 

The  various  expressions  for  Che  component  and  total  forces 
given  In  the  preceding  sections  contain  Che  propeller-appendage  length 
ratio  c  -  c/2i  as  a  parameter.  When  In  these  egressions  we  set  c  »  0, 
the  forces  on  both  propeller  and  appendage  disappear.  However,  In 
spite  of  this  disappearance,  the  dimensionless  ratios  F/T  approach 
definite  non-zero  limits  as  c — ^0.  The  curves  for  the  case  c  »  0 
are  therefore  to  be  Interpreted  as  plots  of  these  limiting  values. 
These  c  -  0  curves  are,.-  .  U  identical  with  the  curves  of 

Breslln  who  obtained  the  forces  on  the  appendage  by  making  the 
simplifying  assumption  that  the  propeller  has  both  a  constant  non¬ 
zero  circulation  and  zero  chord. 

In  Che  discussion  to  follow  we  shall  first  describe  the 
effect  of  varying  x^,  r  and  c  on  the  nature  and  magnitude  of  Che 
hydrodynamic  forces  and  Chen  draw  some  conclusions  with  regard  to 
optimum  design  features  of  propeller-appendage  assemblies. 
feidde^^o8ition(joc^J^ 

Effects  on  Appendage  -  Figures  13  through  16  show  sample  plots 
of  the  component  forces  Indicating  that  while  the  apparent- mass  forces 
go  up  more  than  proportionately  with  c,  the  rr  .  i  i  wake 

forces  change  less  Chan  proportionately  with  c.  The  complex  way  In 
which  these  component  forces  make  up  the  total  force  Is  shown  In  a 
sample  plot  of  i  .  17.  Figs.  19  through  29  show  the  variations  of  the 
total  forces  and  these  can  be  summarized  as  follows: 
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a)  The  forces  are  predominantly  positive,  i  e  ,  they  act 

In  a  direction  opposite  to  propeller  motion.  Fee  reasonable  spacings, 
X  >  1.1,  a  minimum  in  the  value  of  F  /T  occurs  around  c  ---  ,15  for 

0  d 

r  »  2  and  around  c  =  .1  for  r  =  4. 

b)  The  maximum  force  occurs  before  the  propeller  center  passes 

the  appendage,  with  y^^  moving  closer  towards  0  with  a 

decrease  in  c. 


c)  After  the  propeller  has  passed  the  appendage  the  forces 
are  usually  negative  and  negligible 

d)  The  maximum  forces  decrease  with  an  increase  in  propeller 
velocity. is  reduced  by  about  257.  in  going  from  r  -  2  to  r  -  4 

e)  As  shown  in  Figs.  31  and  33  the  forces  go  up  drastically 

with  a  decrease  in  x  .  F^  increases  2-1/2  times  in  going  from 

o  am 

x^  =  1.6  to  x^  =  1.2. 
o  o 

Effects  on  Propeller 

The  apparent  mass  forces  and  also  usually  the  wake  force 
go  up,  according  to  Figs  35  through  37,  more  than  proportionately 
with  c  while  the  quasi  steady  forces  go  up  less  than  proportionately 
with  c.  Fig.  38  shows  how  the  three  component  forces  make  up  the 
total  force.  It  should  be  noticed  that  the  line  Fa/T  ^  1  represents 
the  force  on  the  propeller  at  infinity  or  the  unperturbed  thrust 
Thus  the  plots  in  Figs.  40  to  45  in  addition  to  the  actual  forces  on 
the  propeller  also  represent  the  ratio  of  the  actual  force  to  that 
without  interference.  The  variation  of  the  total  force  can  be  viewed 


as  a  small  or  moderate  perturbation  round  Fp^ /T  = 
propeller  force  always  acts  in  the  same  direction 


1  and  thus  the 
Fp  is  independent 
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of  r  and  its  dependence  on  and  c  can  be  summarized  as  follows: 

a)  The  effect  of  appendage  interference  is  in  general  to  yield 
forces  Liss  than  the  unperturbed  propeller  thrust.  This  reduction 
increases  with  c. 

b)  The  propeller  force  reaches  a  minimum  before  it  passes 
the  appendage  in  the  region  of  0^.y^^,w6 

c)  A  measurable  maximum  as  shown  in  Fig.  47,  is  reached  by 
F  only  at  X  <1.2. 

d)  After  the  propeller  has  passed  the  appendage  the  forces 
on  it  are  slightly  higher  than  the  unperturbed  thrust. 

From  an  analytical  point  of  view  it  is  interesting 

to  note  in  Figures  49  and  51  the  rather  wide  scatter  of  the  coordinates 

of  the  substitution  vortices  The  value  of  x^  even  fails  to  approach 

the  expected  quarter-chord  position  as  y^  — ►  oo 
•  « 

^ke^^ositiot^^-x^^^)^ 

Effects  on  Appendage  -  For  the  case  where  the  propeller  is 

upstream  of  the  appendage  both  the  quasi-steady  forces  and  the  wake 

forces  increase  more  than  proportionately  with  an  increase  in  'c. 

However,  the  apparent  mass  forces  exhibit  a  fairly  complex  behaviour. 

The  depressions  occuring  at  z.  0  were  unexpected  and  their 

influence  can  be  noticed  in  the  shape  of  the  total  force  as  exemplified 

in  Fig.  18.  From  the  succeeding  plots  of  the  total  appendage  force 

for  a  skeg  position  the  following  general  comments  can  be  made. 

a)  Except  for  x^  ^  -1.2  the  forces  are  predominantly  negative 

acting  in  the  direction  of  propeller  motion  The  ratio  F  /T  seems  to 

d 
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1 


Increase  with  an  increase  in  c.  Thus  with  an  increase  in  propellei: 
chord  the  value  of  F.  will  increase  both  by  vin  arc  of  an  increase 
in  c  as  well  as  an  increase  in  the  coefficient  F  /T. 

b)  The  tnaxiouffl  (negative)  force  occurs  before  the  propeller 
center  passes  the  appendage  and  this  occurs  in  the  fairly  narrow 
range,  •!<  Vq  <*5. 

c)  After  the  propeller  has  passed  the  appendage  the  forces 
are  usually  positive  and  negligible. 

d)  The  maxiinunj  (negative)  forces  decrease  with  an  increase 
in  propeller  velocity. 

e)  The  forces  go  up  drastically  with  a  decrease  in  x^. 

f)  The  peak  forces  for  a  skeg  arrangement  are  quantitatively 
about  the  same  and  of  opposite  sign  as  Che  peak  forces  or  a  ' .  1  ler 
arrangement . 

Effects  on  Propeller  -  The  apparent  mass  forces  go  up  less 
than  proportionately  with  an  increase  in  "Z  while  the  quasi-steady 
and  wake  forces  go  up  mor.e  than  linearly  with  Z.  From  Figs.  40  to 
45  the  following  emerges: 

a)  The  effect  of  appendage  interference  is  to  increase  the 
propeller  force  above  that  of  the  unperturbed  thrust. 


b)  As  shown  in  Fig.  48b  F^  reaches  a  maximum  in  the  region 


•3  <  Fq  <  1-0 


c)  As  shown  in  Fig.  47  the  value  of  the  maximum  force 


increases  with  a  decrease  in  x^  and  an  increase  in  c. 
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d)  After  the  propeller  has  passed  the  appendage  the  forces 
on  it  are  slightly  below  the  value  of  the  unperturbed  thrust. 

As  shown  in  Fig.  50  the  scatter  of  the  substitution 
vortices  is  here  smaller  than  for  the  propeller-rudder  arrangement  and 
for  large  values  of  they  approach  the  quarter-chord  position  on 
both  propeller  and  appendage. 
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X.  CONCLUSIONS 

1.  Spacing  closer  than  «  |l . zj between  propeller  and 

appendage  should  be  avoided  for  at  least  two  reasons;  at  |l.2| 

there  are  two  cycles  of  stress  reversals;  and  the  peak  forces  become 
excessive. 

2.  For  a  given  propeller  chord  size  c  there  exists  an  optimum 
ratio  c  which  will  produce  the  smallest  peak  force  on  the  appendage. 
These  values  are  2i  7c  for  low  propeller  speeds  (r  =  2);  and  2i  »10c 
for  high  propeller  speeds  (r  =  4) .  The  above,  of  course,  applies  to 
the  practical  range  of  c  values^  .1<^<.5  and  excludes  the  case  of 

i  s  0  when  the  force  on  the  appendage  would  be  zero. 

3.  Peak  forces  cn  both  appendage  and  propeller  always  occur 
before  the  propeller  center  passes  the  appendage;  the  forces  are 
negligible  after  the  propeller  has  passed  the  appendage. 

4  The  peak  forces  on  the  appendage  do  not  differ  appreciably 
for  either  skeg  or  rudder  arrangement. 

3.  The  peak  forces  on  the  propeller  for  a  rudder  arrangement 
are  essentially  equal  to  the  unperturbed  thrust;  for  a  skeg  arrange¬ 
ment  they  are  measurably  higher. 

6.  On  the  appendage  the  peak  forces  Increase  less  than  pro¬ 
portionately  with  an  Increase  in  V;  on  the  propeller  they  are  pro¬ 
portional  to  Its  velocity. 
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APPENDIX  A 


EVALUATION  OP  INTEGRALS 
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1.  Solution  of 


dn 


(A-1) 


having  a  singularity  at  n  “  ^2-  and  ^2  are  complex  constants. 
Writing 


(n-k2)(n"iii)(n“n2)  n-k2  n-ni  n*n2 


we  have 

A 


kj^-fk2 _ 

(k2“'l2^)  (k2“'n2^ 


B 


C 


..  ■  . - 


Thus  the  first  Integral  of  1^^ 


-±L, 

H-k2 


becomes 


This  Is  a  principal  value  Integral  whose  value  Is  v. 
The  second  and  third  terms  of  are  of  the  form 


(A-2) 
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A2. 


(A- 3) 


with 


p-a+lb  a>i. 

Multiplying  top  and  bottom  of  (A-3)  by  4-  and  writing 


s  -  i  cos  d 


n  -  i+ik 


.  f''  (i+cqselde  .  1  -I  p  .„)de 

cose  +  p  j/jj  cose  +  p 


or 


Il2  "  IT  +  (1  p) 


de 


cose  t  p 


The  problem  Is  now  to  Integrate  the  last  Integral  rewritten  as 


1 

7 


de 

cose  t  p 


(A-4) 


This  Integral  can  best  be  Integrated  In  the  complex  plane  by  a 

le 

contour  Integration  on  the  unit  circle  z  ■  e  .  In  this 
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notation  we  have 

C0.9  -  ^-2^  d9  .  If 

and 

T  _  _ 

^13  ■  tr  j  (z-zj^)  (Z-Z2) 

where 

a)  for  +  p 


with  Zj^  inside  tlie  unit  circle, 
b)  for  “p 


z^  =  p  + 


Z2  - 


P 


(A- 5) 


witii  ^2  inside  tne  unit  circle 
Since  the  values  oi  the  residue 


for  1^3^  is 


2iTi 

Z1-Z2 


and  for 


I 


13- 


2Tri 

Z2-Z1 


we  have 


1 


d9 _ 

cos9  +  p 


(A-6) 
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Also 


In  terms  of  our  parameters 


If  In  addition  and  t)2  are  complex  conjugates,  then  we  have  with¬ 
in  the  brackets  an  expression  of  the  form  ^-^2  ~  ^1^2  since 

^1^2 


h  - 


Ij"'! 


\U  +1 
-^lU  -I 


(A-9) 


2 .  Solution  of 


where 

*2  *  ^o  ■*■ 
ai  -  Yo  -  c/2. 
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Writing 


y-y 


-1 

where 


•ii  ■  I  *12  *  I  tyo-i(v^>' 


(A-ll) 


•13  - 1  lyo-^tv*)’- 

Multiplying  top  and  bottom  by  i/ (1-0  (Ah+C)  »nd  breaking  up  tne 
numerator  into  partial  fractions,  we  obtain  for  I2- 

+  l®ir*13(^"*ll'^*13^^ 

1 

_ ^L— 

^  (ai3+C)p 


-  -  I21  +  ^22  ^*ir*13^^**ll'^*13^^  ^3 


where 


P 


fih  (a^i^+O  (ai2'*-^). 
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We  shall  now  attempt  to  transform  I2  Into  a  function  of  elliptic 
Integrals  by  use  of  the  following  transformation 

.  _  cos^-v 

?  l-VCOljl 


<iC 


(v^-l) 


— Sisi—  d# 

(l-VCOS^) 


where  v  Is  a  constant  to  be  determined  In  such  a  manner  as  to  yield 
for  l2^,  l22»  ^23  Integrals. 

Taking  I22  first,  we  have 


[(v^-l)sln^/(l«vcoe^)^ld^ 


{,  /COS^»V  .  cos^-v  N /.  .  cos^-v  V  I 

“^l-vcos^'  ^®11  ^  1-VCOS0^  ^®12  l-vcos0^ 


T 

7 


l-a',Vv  l-a,«v 

— cos^)(l  +  -  -i'  ■  ■  cos^)] 

®11  ^  h2  ^ 


m 


[(1  + 


The  value  of  v  now  has  to  be  chosen  so  as  to  make  the  coefficient 
of  cos^  In  the  denominator  vanish;  hence 


a^l-v 


+ 


ai2-v 


0 


or 
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By  writing  0  ■  0  +  'TT/l  and  realizing  that  we  are  dealing  with  an  even 
function  the  above  integral  is  recognized  as  twice  the  complete 
elliptic  integral  of  the  1st  kind,  or 
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_ 


J  K(k)  + 


Multiplying  top  and  bottom  of  12^]^  by  (1  +  «)  co80) 

yr  IT 

/»  r 


(1+  i)co8i)di 

J  (1-;>2co8^0)  \/l+k^CO8^0 

o 


cosj  di 


o 


(1-  1>^CO8^0)  v/l- 


k^CO8^0 


The  first  of  the  above  integrals  is  twice  the  complete  elliptic 
integral  of  the  3rd  kind.  The  second  integral  being  odd  in  0  is 
zero .  Thus 


Considering  next  I23 


'rr 

/* 


'■23 


1-^ 


>> )  (ai2"  ^  ^ 


di 


o 


1*13^  cosi- »)  1  |/l-k‘co8  0 


(A-13) 
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vCa, --v) 

+  [1  +  _i3 - ] 

^  i«va  * 


13 


[1 


_ a . .  I 

+  -  -  rP  cos^J  \/l-k^co8^^  J 
®13  ^ 


Multiplying  the  second  Integral  top  and  bottom  by 
l-va. 


1  - 


‘13 


ai3-v 


COS^ 


L23 


— 


f  2v(a,--v) 


*13“’' 


A9. 


+ 

T^va 


13 


l“vai  <j 


*13“’' 


o  [1-  (  cos^^]  \J l-k^coB^0 


] 


or  since  the  second  part  of  the  second  Integral  Is  zero 

-7 - 


^3  “  T7a 


13 


l*y 

(*ll-v)(*,,-v) 


2  ‘ 

with  -  (-r— r~) 
*13  ’' 


The  total  expression  for  I2  Is  then 

T  -  2a+v)  /  1-v^ 

2  (l-va^3)  \^l  (*ii-v)(a][2-v) 

l^(l-a^jv)  K(k)  +  (l-v)(l-va^3)  7f  (v^.k) 

v(l”v)  (*^  ^“*^3)  (l”ai~)  if  2  ? 

+ - ^  I'  >w; 


(A- 14) 


(A- 15) 


This  being  a  solution  of  Integral  (A- 10)  In  terms  of  complete  elliptic 
Integrals  of  the  1st  and  3rd  kind. 
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APPENDIX  B 

SUBSTITUTION  VORTEX  FOR  AEROFOIL  IN  FIEIJ)  OF  ANOTHER  VORTEX 
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Bl. 


Figure  B.l 

Let  there  be  a  vortex  of  strength  (circulation)  Tp 
at  Zp.  Suppose  the  flow  associated  with  the  vortex  in  the  pre¬ 
sence  of  the  aerofoil  (a  flat  plate)  has  circulation  around 
the  latter.  We  seek  the  location,  z.,  of  that  vortex  by  which 
the  aerofoil  may  be  replaced,  such  that  the  flow  at  sufficiently 
large  z  due  to  the  two  vortices  approximates  the  original  flow. 
More  precisely,  z^  is  to  be  so  determined  that  the  flow  due  to 
the  two  vortices  agreijs  with  the  original  flow  to  terms  of  order 

1 

z‘ 

The  configuration  in  question  is  shown  in  Figure  B.l, 
where  it  will  be  seen  that  we  have  chosen  the  appendage  as  the 
aerofoil  which  is  to  be  replaced  by  a  vortex.  By  a  treatment 
entirely  analogous  to  the  one  we  shall  describe,  the  propeller 
may  be  similarly  replaced,  but  it  will  suffice  to  explain  the 
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B2. 

method  with  reference  to  the  appendage  and  merely  quote  the  results 
for  the  propeller. 

Let  w^  be  the  complex  potential  of  the  flow  associated 
with  the  two  vortices.  Then 

"l  ■  W  ?a 

Expanding  the  second  term  In  Inverse  powers  of  z,  we 

have 

''l  "  ^a  ^  ta  -T  +  “^z^ 

where  the  remainder  Is  a  power  series  In  ^  convergent  for 

1*1  >I*J  . 

Now  denote  by  W2  the  potential  of  the  actual  flow 
Induced  by  the  vortex  located  at  in  the  presence  of  the  aero¬ 
foil.  We  first  express  W2  as  a  function  of  where  C  is  defined 

by  the  transformation 

,2 

B.3)  z  -  C  +  ^ 

In  the  H,  plane  the  aerofoil  is  mapped  into  a  circle 
of  radius  with  center  at  0.  The  vortex  at  maps  into  a 
vortex  at  with  circulation  Likewise  the  circulation 

around  the  aerofoil  is  preserved. 

(7) 

Using  Mllne-Thompson' s  circle  theorem  in  the  ^  plane, 

we  get 

''2  "  ^  log  -  Tp)  +  ^f-'^iog  C 
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where  f  Is  the  Image  of  and  the  bar  over  a  letter  stands  for 
P  P 

the  complex  conjugate. 

We  now  expresa  W2  as  a  function  of  z  In  the  following 

formt 

(B.5)  W2  -  +  f(*) 

where  f(z)  Is  a  function  analytic  outside  the  circle  z  ■  i. 

We  then  expand  f(z)  In  Inverse  powers  of  z  as  far  as  the  ^  term. 

z 

To  this  end,  we  Invert  the  relationship  (B.3)  to  obtain 
(B.6)  ^  -  I  +  ^ 

(B. 7)  Cp  -  +  -r 

where  the  sign  of  the  square  root  Is  such  that  as  |z)-^  ». 

We  then  have 
<B.8)  <  -  I  -  +  t.  1) 

(B.9) 

The  remainder  terms  In  (B.8)  and  (B.9)  are  power  series 

In  -7  which  converge  for  (2.|>i. 
z 

We  also  have 

(B.io)  ^  ^  1 

It  Is  readily  shown  that  the  right  side  of  (B.IO)  can 
be  expanded  to  give  the  following  result: 
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N|)-> 


(B.ll) 


B4. 


C-Cp  -  (*-*p)(l  +  ^  7  + 


Substituting  from  Situations  (B.8)  through  (B.ll) 
Into  (B.4),  we  get 


in  ^ 

(B.12)  Wj  -  log  [(z-Ip)(l  +  i  +  o(i))  ) 


-^iogr-Tpa-^i>o(i))i 

ir  a  *2  , 

+  iog  [z(l  -  — ^  o(— t))  ] 

4z^  tT 


We  expand  the  right  side  of  (B.12)  as  far  as  terms 
In  ^  and  obtain,  after  neglecting  an  Irrelevant  constant 

(B.13)  wj  -  >2^  fog(z-Zp)  +  z 

where  the  remainder  Is  a  power  series  convergent  for  |zi>  1. 

Comparison  of  (B.13)  with  (B.2)  to  terms  of  order 

now  gives 


from  which  follows 
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(B.15) 


Equation  (B.IS)  expresses  the  location  of  the  vortex 

replacing  the  appendage  In  terms  of  (since  t  Is  related  to  z. 

r  P  P  P 

by  (B.7)),  and  the  ratio  of  the  propeller  circulation  to  the 

‘oa 

circulation  around  the  appendage.  It  Is  Important  to  note  that 
since  Tp  and  are  real,  (B.15)  shows  that  z^  Is  real  and  that 
therefore  the  substitution  vortex  for  the  appendage  lies  on  the 
appendage  regardless  of  the  location  z^  of  the  other  vortex.  It 
Is  likewise  true  that  the  substitution  vortex  for  the  propeller 
lies  on  the  propeller. 

r 

To  determine  In  (B.15),  we  apply  the  Kutta  condition 
*oa 

at  the  trailing  edge  of  the  appendage,  l.e.,  at  z  -  I.  The  trail¬ 
ing  edge  maps  Into  the  point  ^  In  the  plane  so  that  the  Kutta 
condition  can  be  written  as 

dwo  - 

(B.16)  -  0,  C  “ 


where  W2  Is  given  by  (B.4) 


(B.17) 


Carrying  out  the  di f fprontiation  we  obtain 
dw«  ir„_  f  .2 


^2  _  "‘op,  1  1  . 


in 


oa  1 
w  z 


Application  Of  (B.16)  yields 


(B.18) 


OP-  .  2 
Toa  i 
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Which  when  coiid>lned  with  (B.15)  gives 


(B.19) 


'a  "  i 


1  4.  1 

rz  rr. 


Equation  (B.19)  la  the  desired  explicit  relationship 

between  the  location  of  the  appendage  vortex,  and  z^,  the 

location  of  an  arbitrary  "exciting"  vortex.  By  way  of  checking 

this  result,  we  observe  that  as  z^  (and  therefore  Cp)  *  ■*> 

z^  ♦  -  ^  ,  the  point  on  the  appendage  distant  ^  chord  from  the 

leading  edge.  This  Is  to  be  expected  since  the  Incident  flow 

approaches  a  uniform  stream  as  z^  »^and  for  a  uniform  incident 

stream  at  any  angle  of  attack  other  than  zero,  the  substitution 

vortex  replacing  a  thin  aerofoil  is  known^^^  to  be  located  at 
1 

the  2^  chord  point.  An  exception  arises  if  »  in  a  direction 
parallel  to  the  y  axis,  in  which  case  the  substitution  vortex 
position  depends  on  x^,  the  distance  from  the  y  axis  to  the 
straight  line  along  which  z^  -»■  <x>.  The  reason  for  this  is  that 
two  limiting  operations  are  not  Interchangeable,  viz  ,  allowing 
Zp  to  00  in  an  arbitrary  direction  on  the  one  hand,  and  allowing 
that  direction  to  approach  perpendicularity  to  the  appendage  on 
the  other.  If  the  latter  operation  is  performed  first,  the 
result  might  be  expected  to  be  anomalous  since  in  this  case  the 
limiting  uniform  stream  Impinges  on  the  appendage  with  a  zero 
angle  of  attack  and  therefore  Induces  no  circulation. 
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To  dotermln*  th«  position  of  the  substitution  vortex 
for  the  propeller,  we  proceed  In  s  manner  exactly  analogous  to 
the  above,  except  that  we  must  remeinber  that  the  Incident  flow 
on  the  propeller  comes  not  only  from  a  vortex  at  z.  but  also  from 
a  uniform  stream  traveling  with  velocity  U  In  the  positive  x  direc¬ 
tion.  With  this  minor  modification  we  get  the  following  expression 
for  Zp,  the  location  of  the  substitution  vortex  that  replaces 
the  propeller: 


(B.20) 


where  *o  ■  +  ly^i  tbe  position  of  the  center  of  the  propeller, 

c  Is  the  propeller  chord,  and  Cl  Is  defined  by 


(B.21) 


Thus  Equation  (B.20)  Is  the  analogue  of  Equation  (B.15) 

In  that  It  expresses  In  terms  of  z^,  r,  :Q^»and  r^p*  Note  that 
since  the  last  two  terms  on  the  rl^t  of  (B.20)  are  Imaginary, 

Zp  lies  on  the  propeller. 

By  application  of  the  Kutta  condition  at  the  trailing 

edge  of  the  propeller,  we  could  derive  an  equation  analogous  to 

(B.ld).  Such  an  equation,  combined  with  Equations  (B.15),  (B.18) 

(or  B.19),  and  (B.20)  would  give  four  equations  for  the  determination 

of  the  four  unknowns,  end  z_.  However,  we  already 

op  os  p  a 
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have  a  fourth  equation,  viz.,  Equation  (6*8)  for  r^ln  the  main 
text  of  the  report,  this  having  been  derived  In  an  Independent 
way  which  also  Invokes  the  Kutta  condition  at  the  trailing  edge 
of  the  propeller.  Thus,  no  further  relationships  are  required. 

One  final  remark  should  be  made  concerning  a  llmlta* 
tlon  on  the  use  of  the  substitution  vortex  method,  To  terms  of 
order  the  potential  due  to  the  presence  of  the  appendage  is 
equal  to  that  due  to  the  appendage  vertex,  but  this  equality 
holds  only  for  observation  points  z  for  which  le,,  obser¬ 

vation  points  exterior  to  a  circle  centered  on  the  appendage  and 
having  a  radius  equal  to  the  half  chord  Similarly,  the  potential 
due  to  the  propeller  is  approximated  by  that  of  the  propeller 
vortex  for  points  outside  a  corresponding  circle  centered  on  the 
propeller.  It  follows  that  In  creating  the  flow  with  both  aero¬ 
foils  present  by  the  substitution  vortex  method,  neither  aerofoil 
should  penetrate  the  circle  associated  with  the  other  This  re¬ 
quirement  Imposes  a  restriction  on  the  variables  t,  c,  and  x^,  viz. 

(B.22)  Xo  >  <  +  f  or  -if.  \J-y\  .  (Xq-I)^ 

Thus  the  method  breaks  down  if  tne  propeller  passes 
too  close  to  the  end  of  the  appendage  On  the  other  hand,  it 
seems  likely  that  a  slight  penetration  of  the  circles  by  the  aero¬ 
foils  will  not  seriously  spoil  the  accuracy  of  the  method  so  that 
practically  speaking,  the  Inequality  (B  22)  Is  a  little  too 
restrictive. 
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